Conformal anomaly as a source of soft photons in heavy ion collisions 
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We introduce a novel photon production mechanism stemming from the conformal anomaly of 
QCDxQED and the existence of strong (electro) magnetic fields in heavy ion collisions. Using 
the hydrodynamical description of the bulk modes of QCD plasma, we show that this mechanism 
leads to the photon production yield that is comparable to the yield from conventional sources. 
This mechanism also provides a significant positive contribution to the azimuthal anisotropy of 
photons, 112, as well as to the radial "flow". We compare our results to the data from the PHENIX 
Collaboration. 
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The manifestations of quantum anomalies in the collec- 
tive dynamics of relativistic plasmas have attracted con- 
siderable interest recently. In the presence of background 
fields, triangle anomalies can lead to non-conservation of 
currents. In the case of the axial AVV anomaly involv- 
ing an axial and two vector currents, the presence 
of an external magnetic field and a finite density of chi- 
ral charge leads to the generation of electric current in 
QCD plasma - the chiral magnetic effect ITWSl. At finite 
density of (vector) baryon charge, magnetic field induces 
the flow of axial current; this is the chiral separation ef- 
fect [SHI]. Both of these phenomena are an integral part 
of relativistic hydrodynamics, and in fact are required by 
the second law of thermodynamics [5HT1J . 

In this letter we investigate the related effects stem- 
ming from the conformal SVV anomaly [121 that in- 
volves a scale (dilatational) current S'p and two vector 
currents and reflects the violation of conformal in- 
variance of QCD by quantum effects. The conformal 
anomaly results from the running of the coupling con- 
stant (asymptotic freedom [13 ) in QCD and expresses 
the non-conservation of the dilatational current S'^, so 
that the trace of the energy-momentum tensor 0^ does 
not vanish even in the chiral limit of massless quarks: 
d'^Sfi — 9^. The quarks carry both color and electric 
charges, so when QCD is coupled to electromagnetism, 
the quark triangle diagram induces an anomalous cou- 
pling of the trace of the energy-momentum tensor to 
photons |14H16j . The trace of the energy-momentum 
tensor in hydrodynamics excites the bulk modes of the 
fluid that are abundant in (non-conformal) quark-gluon 
plasma [T7lf20] . The heavy ion collisions at early times 
produce very strong background magnetic fields O |2T] . 
As a result, the conformal anomaly acts as a source of 
photon production that is powered by the energy of the 
bulk hydrodynamical modes in the plasma. This is the 
mechanism of photon production that will be discussed 
in detail below. Note that while we will use hydrody- 
namics to describe the bulk modes in the plasma, the 
deviation from equilibrium in general need not be small 



for our mechanism to operate. For example, the non- 
equilibrated Bose-Einstein condensate of gluons [22] may 
be even more effective in producing photons. Note that 
unlike in the conventional scenario, the quarks in our case 
appear only in the triangle loop that receives contribu- 
tions from the virtual UV modes - so the production of 
real on-shcU quarks is not required, and the mechanism 
can operate even at very early times. 

Because of the relatively weak interactions with the 
medium, "direct" (i.e. not resulting from the hadron 
decays) photons play an important role of a "thermome- 
ter" of the quark-gluon plasma [23j [24] since the rate 
of their production per unit volume is expected to scale 
with the temperature T as ^ T^. Recent measurements 
of the PHENIX collaboration show very large excess of 
direct photons with the transverse momentum up to 3 
GeV in AuAu collisions at RHIC [25]. The azimuthal 
anisotropy (often called "elliptic flow") of the produced 
photons has also been measured and reported [25^. It 
has been found that the anisotropy is large and similar 
to the elliptic flow of hadrons. This result contradicts the 
current theory of photon production. Indeed, it has been 
expected (see e.g. [26]) that the elliptic flow of direct 
photons would be much smaller than that of hadrons be- 
cause a significant fraction of photons has to be produced 
at early times, when the temperature is the highest. The 
measured yield of soft photons is indeed large and re- 
quires an early time production mechanism (see [271 128j 
and references therein). At these early times the hydro- 
dynamical flow has not been built up yet, and so the 
photons produced at that time are not expected to pos- 
sess a significant azimuthal asymmetry [57]. Moreover, 
the jet-medium interactions and the resulting induced 
bremsstrahlung of photons is expected to lead to a neg- 
ative contribution to the elliptic flow Fourier coefficient 
V2 [22] - indeed, due to the geometry of the collision the 
produced medium has an almond shape and is elongated 
along the axis orthogonal to the reaction plane. 

Therefore the large and positive value of V2 presents 
a serious challenge to theory. Here we will consider the 
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FIG. 1: The coupling of the conformal anomaly to the ex- 
ternal magnetic field resulting in photon production. 



photon production mechanism stemming from the con- 
formal anomaly of QCD x QED and the presence of a high 
magnetic field in heavy ion collisions. We will demon- 
strate that this mechanism results in a significant pho- 
ton and dilepton yields that are comparable to the ones 
from the "conventional" mechanism and may potentially 
explain the V2 puzzle for soft direct photons. 

Let us begin by reminding the basics of conformal 
anomaly. In field theory the divergence of the dilatational 
current is equal to the trace of the energy-momentum 
tensor. In QCD, this divergence does not vanish signal- 
ing the breaking of scale invariance due to dimensional 
transmutation and the running coupling: 
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G^''"G^.a+^™, [1+ 7m (5)] 99, 



(1) 

where /3{g) is the beta- function of QCD, rriq are the quark 
masses, and 7m (ff) Eire the corresponding anomalous di- 
mensions. The current acting on the vacuum produces 
scalar color-singlet states a of mass TOq. with an ampli- 
tude fa'. 



{0\S^\a)=iq'^U; {0\d,S'^\a)=ml 



(2) 



Let us now consider the coupling of QCD scale 
anomaly to electromagnetism. This coupling can lead 
to the production of photons in external magnetic field 
as described by the diagram of Fig. [T] To evaluate the 
contribution of this diagram, we need to consider the 
coupling of the scalar meson to photons. This coupling 
is described by the triangle quark diagram, and leads to 
the following effective interaction 14-16J: 
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where 50-77 is related to the decay constant /o- discussed 
above and to the ratio of cross sections of e^e~ annihi- 



lation into hadrons and muons 
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where a is the fine structure constant. The resulting 
width of a decay into two photons is given by P^IHTB] 



r((T 77) = gl 
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Using R — 5 for six quark flavors (all of which contribute 
to the triangle diagram) and the values nia- = 550 MeV, 
fa = 100 MeV discussed above, we get from ^ the value 
r(CT 77) ~ 5 KeV. This is in the middle of the range 
(2 10 KeV) for the two photon decay width of /o(600) 
meson listed by PDG [SOj, supporting the identification 
of the lightest a dilaton with this meson. This allows us 
to fix the value 50-77 — 0.02 GeV~^. Now we have all the 
information necessary to evaluate the diagram of Fig. [T] 

To compute the photon production rate from the di- 
agram of Fig. [1] we evaluate the imaginary part for the 
photon self-energy, see [SI] [32] . A straightforward cal- 
culation yields for the production rate at mid-rapidity 
(Qz = 0) the following expression: 
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peiqo = \q\). (7) 



Since we consider production of photons in the QCD 
plasma, it is appropriate to use the hydrodynamic spec- 
tral function of the bulk mode 9 [331 [311 : 
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where Tg = (4/3?7 -f- C)/(e + p) is the sound attenuation 
length and r/ and ( are shear and bulk viscosities. The 
second term describes the sound peak at qo = Cs\q\. The 
sound mode does not contribute to the production of real 
photons since the width of the sound peak is not large 
enough to reach the null dispersion of photons. Therefore 
the photon production is dominated by the bulk viscosity 
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In deriving Eq. ([t]) we neglected the z-component of 
the magnetic field, because it is expected to be an order 
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of magnitude smaller than and By {B^ ~ Bx.yl^)\ we 
also neglect the contribution of the electric field. 

In what follows we will compare our result with the 
baseline provided by the conventional thermal photon 
production rate recently calculated on lattice [35]: 
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where Cem = ij"^ = J2fQf with Q/'s are the electric 
charges of the quarks, and pv is the vector current spec- 
tral function that in the limit of go and q — > is 
related to the electric conductivity: 
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Note that this conventional mechanism ( 10 ) is expected 



to be the dominant one for low transverse momentum, 
p±, photons. For photons with p± ^ 2 GeV and above 
there will be additional contributions to the rate which 
can be calculated perturbatively. However we did not 
include these additional contributions as we are mainly 
interested in low p± photons. 

The spectral function for 9 and the bulk viscosity was 
calculated in lattice QCD [El |33] . However the extrac- 
tion of bulk viscosity from the lattice data is notoriously 
difficult. To get an independent estimate of the bulk vis- 
cosity we thus follow |3ni 133 and assume that 

2 



c 



(12) 



Thus the bulk viscosity vanishes in the conformal limit, 
= 1/3. In the relaxation time approximation, this ex- 
pression is obtained in the kinetic theory with = 15 
(see e.g. [40]). The paper [40] contains also a phenomeno- 
logical estimate of the value of bulk viscosity inferred 
from the comparison of viscous hydrodynamical compu- 
tations with the data on the elliptic flow of mesons and 
baryons. The resulting estimate is (/s = 0.005 [IQ. Us- 
ing the lattice data for the speed of sound in the freeze- 



out temperature range from Ref. 



0.175H-0.221, 



we infer for the bulk viscosity from ( 12 ) the value of 



~ 2.5 5. The leading log calculations in SU(3) Yang 
Mills theory results in a much larger value ~ 48, see 
Ref. [30] • In our calculations, we choose the lowest value 
available in the literature, — 2.5-;- 5, with an assump- 
tion rj/s = l/47r. 

The magnetic field in heavy ion collisions was esti- 
mated in Refs. [3] and [21]; the fluctuations of magnetic 
fleld were evaluated in Refs. [38] and [39]. In this paper, 
we neglect the spatial gradients of magnetic field and 
estimate the time dependence in the eikonal approxima- 
tion taking into account only the (leading at large times) 
contribution from spectators: 
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FIG. 2: The azimuthal anisotropy V2 of the direct photons 
for different values of bulk viscosity corresponding to Cq in the 
range of 2.5-^5 calculated for minimum bias Au-Au collisions. 
The dashed line represents the results with = 4. The black 
dotes are the data from the PHENIX collaboration [25] for 
minimum bias Au-Au collisions at y/s — 200 GeV. 



where eBf it the magnitude of the i-th component of the 
magnetic field a,t t — and tg is the characteristic decay 
time. The a;-component of magnetic field at i = 0, B^, 
is approximately independent of the impact parameter 6, 
while the y-component is linear in b. Both components 
B'^ y are linear as a function of the collision energy, ^/s^, 
the typical decay time is inversely proportional to ^/s. 

Here we neglect the transverse expansion of the fire- 
ball and assume that it has an almond shape with the 
following characteristic sizes in x and y direction: — 



{RA — b/2) and ly = \/ R\ — where Ra is the radius 
of the colliding nuclei. We approximate the time evolu- 
tion of the temperature at early times using the Bjorken 
hydrodynamics T/Tq = (tq/t)^/^, where Tq is the initial 
temperature and tq is the initial time (given by the char- 
acteristic thermalization time of the gluons) that can be 
estimated in terms of the saturation scale, Qs, and the 
coupling constant, as, see e.g. Ref. [22]. For Au-Au col- 
lisions at y/s = 200 GeV we use tq = 0.1 fm/c. 



To evaluate the bulk viscosity (12) we need the speed 
of sound, Cs and the entropy, s; we use the model 
parametrization ,42^ of lattice results for pure glue SU(3) 
theory. Note that the transport coefficients of the plasma 
may be affected by magnetic field; for recent examples, 
see (43] and [44] . 

Our results for the azimuthal anisotropy of photons 
calculated using both conventional production mecha- 
nism and the one from the conformal anomaly are shown 
in Fig. [2] for the minimum bias Au-Au collisions at 
\/s = 200 GeV. In our approximation (no transverse 
fiow) , the conventional mechanism does not give any con- 
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FIG. 3: The transverse momentum spectra of the produced 
direct photons for Cq = 2.5 calculated for minimum bias Au- 
Au coUsisions, see text for details. 



tribution to the azimuthal anisotropy. The comparison 
with the experimental data from PHENIX [55] indicates 
that conformal anomaly could account for a large fraction 
of the observed photon anisotropy. 

In Fig. [3] we show our result for the transverse mo- 
mentum spectrum of direct photons. Due to the factor 
of in the production rate ([7|, the spectrum of pho- 
tons produced due to conformal anomaly is enhanced in 
comparison to the conventional one at transverse mo- 
menta k±> \ GeV. The factor of in the rate hardens 
the transverse momentum spectrum, and magnetic field 
grows with the impact parameter of the collision; these 
two effects thus conspire in mimicking both the elliptic 
and radial flow of photons in non-central collisions. 

An interesting corollary of our mechanism is the po- 
larization of the produced photons (and low-mass dilep- 
ton pairs) relative to the reaction plane of the collision. 
Other tests include the study of U-U collisions, where 
the deformed shape of the U nucleus may allow to sepa- 
rate |45j the eccentricity of the initial condition from the 
magnitude of magnetic field that drives our effect. 

The calculations performed in the current paper are 
quite schematic and rely on a number of crude approxi- 
mations. The bulk viscosity and its temperature depen- 
dence in SU(3) Yang-Mills theory are the major sources 
of uncertainty in our calculations. Nevertheless, in this 
letter we preferred to err on the side of caution and used 
the most conservative estimates for the bulk viscosity and 
other input parameters. In spite of this, we find that the 
quantum anomaly is responsible for a very substantial 
contribution to the overall soft photon yield. Realistic 
calculations treating the 3D hydrodynamical expansion 
and proper initial conditions are required to reach a defi- 
nite conclusion, and to compare to the data on the trans- 



verse momentum spectra of photons. These calculations 
are in progress and will be presented elsewhere. 
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